Raman study of Fano interference in p-type doped silicon 
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As the silicon industry continues to push the limits of device dimensions, tools such as Raman 
spectroscopy are ideal to analyze and characterize the doped silicon channels. The effect of inter- 
valence band transitions on the zone center optical phonon in heavily p-type doped silicon is studied 
by Raman spectroscopy for a wide range of excitation wavelengths extending from the red (632.8 nm) 
into the ultra-violet (325 nm). The asymmetry in the one-phonon Raman lineshape is attributed 
to a Fano interference involving the overlap of a continuum of electronic excitations with a discrete 
phonon state. We identify a transition above and below the one- dimensional critical point (Er x = 3.4 
eV) in the electronic excitation spectrum of silicon. The relationship between the anisotropic silicon 
band structure and the penetration depth is discussed in the context of possible device applications. 
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I. INTRODUCTION 

As the silicon industry continues to push the limits 
of device dimensions, tools such as Raman spectroscopy 
are ideal to analyze and characterize the doped silicon 
channels being developed. Due to the small penetration 
depth (~ 8 nm) of 325 nm light in silicon and motivated 
by recent doped silicon-on-insulator (SOI) devices, con- 
focal micro-Raman spectroscopy offers a unique method 
to probe both bulk and surface scattering. Raman scat- 
tering in solids can arise through a number of different 
processes including scattering by phonons, plasmons, in- 
traband, and interband electronic transitions. Intraband 
electronic Raman scattering, where the initial and final 
states are in the same band and the intermediate state 
is in a different band, has been thoroughly investigated 
both experimentally and theoretically [1]. Only a few 
experiments have been done on interband electronic Ra- 
man scattering and research has focused on inter- valence 
band scattering in p-type doped silicon Q, H, 0, • 

Experiments were conducted to investigate the effect of 
free holes on the optical phonons of heavily p-type doped 
silicon as well as an asymmetry in the Raman lineshape of 
the zone center optical phonon (521 cm -1 ) 0, Q. When 
exposed to incoming light, the doped silicon produces a 
continuum of inter- valence band electronic excitations, 
which overlaps the Raman-active one-phonon energy. As 
a result, the overlap between the electronic continuum 
and discrete state causes interference effects and a phe- 
nomenon termed Fano interference arises. Fano interfer- 
ence was first described by U. Fano after observing an 
asymmetry in the electron energy loss spectrum of he- 
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Hum [6j]. The asymmetry was attributed to a quantum 
interference involving the overlap of the wavefunctions of 
discrete states and a continuum. The helium lineshape 
was fitted with a Fano function, which depends on two 
parameters. In this paper, we present a detailed experi- 
mental and theoretical study of the Raman spectrum of 
p-type doped silicon for a range of carrier concentrations 
and a wider array of excitation wavelengths, extending 
into the ultra-violet (UV), to probe the surface and bulk 
silicon states. 



II. THEORY 

The Fano effect in p-type silicon is associated with the 
resonant interaction of a discrete state with a continuum 
(Fig. 1). The first state, being a discrete state, represents 
a hole in its acceptor ground state and one locally excited 
optical phonon. The second state, being a continuum, 
represents no phonons being created and a hole excited 
into the p 3 / 2 valence band Q. The Fano effect in p- 
type silicon has been observed in photoconductivity and 
Raman spectroscopy experiments 0] • The asymmetry 
due to the zone center (k = 0) optical phonon changes 
as a function of temperature, excitation wavelength, and 
doping concentration [5] . Interaction solely with the hole 
continuum is due to the fact that shallow acceptor levels 
are separated by energies comparable to phonons. 

The changes in the Raman lineshape are due to the ex- 
istence of a continuum of electronic excitations produced 
by inter- valence band transitions whose energy overlaps 
that of the zone center optical phonon. As the boron 
doping concentration is increased in silicon, the light and 
heavy hole bands are split further from one another. One- 
electron excitations, which are infrared forbidden but Ra- 
man allowed at k = 0, scatter electromagnetic radiation 
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FIG. 1: (a) Isotropic parabolic energy bands with the Fermi 
energy Ei? in the valence bands [10]. Diagram of Fano in- 
terference: (b) Discrete state; represents a hole in its accep- 
tor ground state and one locally excited optical phonon (c) 
Continuum; represents no phonons being created and a hole 
excited into the p 3 / 2 valence band. 



which is observed as Raman scattering Q . Using a mod- 
ified Fano function, 
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where 



, Cardona et at. were able to fit the 



experimental Raman data of p-type doped silicon and 
extract the two parameters q and r [f|. T is the squared 
matrix element of the coupling between the continuum 
and the discrete phonon state, which depends only on 
the carrier concentration, while q squared represents the 
ratio of the scattering probability of the discrete state to 
that of the continuum, and depends on both the carrier 
concentration and excitation wavelength [J]. 

The peak position of the intrinsic one-phonon Raman 
peak is shifted as the doping concentration is increased. 
The doped silicon critical frequency (ft) is related to the 
phonon frequency by [H[ 
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where Qi n t represents the intrinsic silicon peaks and 5Q 
represents the shift of the critical frequency from the in- 
trinsic frequency. The shift is related to the parameters 
q and T by @ 



(3) 



where SQ max represents the shift in the position of the 
maximum peak of the doped silicon peak from the intrin- 
sic silicon peak. 

Crystalline semiconductors exhibit first-order Raman 
scattering only by phonons with k = [8] . Since the pho- 
ton wave vector (k p hoton = 2tt/X) is much less than the 
Brillouin-zone size (27r/ao), one can assume vertical tran- 
sitions and set q = and satisfy conservation of crystal 
momentum (k s = ki ± q) [9|. Incoming light of a par- 
ticular wavelength (ui) excites electrons in the light hole 
band to the conduction band, which then interact with a 
phonon, emit a scattered photon (a;') and de-excite into 



the heavy hole band. Assuming isotropic parabolic en- 
ergy bands, these inter- valence band transitions give a 
Raman spectrum with cutoffs at Q, m and Om, represent- 
ing the minimum and maximum energies defined by the 
valence band curvature and the Fermi level for a given 
carrier concentration (Fig. 1) [10]. However, silicon has 
highly anisotropic valence bands and in order to obtain 
quantitative results it is necessary to have knowledge of 
the detailed band structure over the whole Brillouin zone. 

The Fermi level is not dependent on the excitation 
wavelength. As the excitation energy is varied, corre- 
sponding to different energy gaps for the inter-valence 
band transitions, different orientations in &-space will 
have higher Raman efficiencies. Doping may also shift 
the valence and conduction bands vertically as well as 
the Fermi level. For each excitation wavelength, the 3D 
E — k plot for silicon should be analyzed to determine 
which crystal orientation will dominate and which light- 
heavy valence band transitions will contribute to the Ra- 
man spectrum. For certain crystal orientations, such as 
[110], the minimum energy Q m is no longer less than the 
zone center optical phonon (O m > Q = 65 meV), which 
disallows the overlap of the continuum with the discrete 
phonon state [Bj]. 



III. EXPERIMENT AND RESULTS 

Raman data was collected from a triple-axis spectrom- 
eter (Jobin Yvon Horiba, T64000) using a liquid nitro- 
gen cooled CCD detector in the back-scattering config- 
uration. Measurements were performed with several cw 
lines of a HeNe (632.8 nm), diode (532 nm), and HeCd 
(442 nm and 325 nm) laser and the second harmonic of a 
pulsed Tiisapphire laser (Mira 900, Coherent, 5 ps pulses, 
76 MHz repetition rate, 382.6 nm). The excitation light 
was focused onto a sample using a long distance objective 
(Nikon, 50x/0.45) to a spot size of ~ 1 /im. The carrier 
concentrations for the silicon samples were determined 
by sheet resistance and Hall Probe measurements. 



A. Si(100) Samples 

Measurements were performed at room temperature 
on thermally diffused boron doped Si(100) wafers, which 
were cleaned with a buffered oxide etch (BOE) and dried 
with N2. The Raman data obtained for the three ther- 
mally diffused boron doped samples are shown in Figures 
2 and 3. The experimental data (discrete points) were 
then fitted with the modified Fano function (solid lines) 
with a least-squares computer fit which used q and T as 
adjustable parameters. 

The values extracted from the fit for T, <$Q, and 
fi^max are listed in Tables I and II. The fit obtained for 
Fig. 2 has excellent agreement (< 2% error), whereas the 
fitting function begins to deviate from the experimental 
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FIG. 2: Raman lineshapes vs. scattering wavelength for three 
thermally diffused boron doped Si(100) samples at excitation 
wavelengths 632.8 nm, 532 nm, and 442 nm (1800 gr/mm). 
The experimental data (discrete points) were then fitted with 
the modified Fano function (solid lines) with a least-squares 
computer fit which used q and T as adjustable parameters. 
The spectra were normalized by the maximum peak value of 
each concentration. 



values as the excitation wavelength moves deeper into the 
UV. 

As the excitation wavelength decreases towards the 
UV, the T values are no longer comparable and the q 
values become increasingly higher and then negative. If 
T is held constant, which according to Cardona et at. 
should be independent of the excitation wavelength, for 
a doped silicon sample with carrier concentration 2.50 x 
10 20 cm -3 , the experimental data can be fitted as a func- 
tion of q only and compared to the different excitation 
wavelengths (Fig. 4). 
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FIG. 3: Raman lineshapes vs. scattering wavelength for three 
thermally diffused boron doped Si(100) samples at excitation 
wavelengths 382.6 nm and 325 nm (2400 gr/mm). The ex- 
perimental data (discrete points) were then fitted with the 
modified Fano function (solid lines) with a least-squares com- 
puter fit which used q and Y as adjustable parameters. The 
spectra were normalized by the maximum peak value of each 
concentration. 



TABLE I: Parameters obtained by fitting experimental first- 
order Raman spectra with the modified Fano function (pQ). 
Spectra were obtained for three different excitation wave- 
lengths (632.8 nm, 532 nm, and 442 nm). 



Doped Si (cm^; 
Concentration 


T (cm" 1 ) 
632.8 532 442 632.8 


q 

532 


442 


2.50 x 10 2U 

4.97 x 10 19 

6.98 x 10 18 


9.22 9.39 9.12 1.39 
4.05 3.79 3.33 7.71 
1.93 1.98 1.87 23.6 


3.13 8.01 
14.6 35.4 
37.3 32.2 




Doped Si (cm -5 ) 
Concentration 


(cm x ) SQmax (cm 
632.8 532 442 632.8 532 


442 


2.50 x 10^ u 

4.97 x 10 19 

6.98 x 10 18 


-7.8 -7.8 -5.8 -1.17 
-0.8 -2.0 -0.3 -0.27 ■ 
-0.3 -1.6 -0.6 -0.22 ■ 


-4.8 
-1.74 
-1.55 


-4.66 
-0.21 
-0.54 



The values extracted from the fit for g, 5Q, and SQ max 
are listed in Table III. The theoretical fit has excellent 
agreement for all values down to the 382.6 nm excitation, 
but does not converge for the 325 nm excitation. Also, 
the q values are negative for the 325 nm excitation. The 
interpretation of these results will be addressed below. 



TABLE II: Parameters obtained by fitting experimental first- 
order Raman spectra with the modified Fano function JT]). 
Spectra were obtained for two different excitation wavelengths 
(382.6 nm and 325 nm). 



Doped Si (cm" a ) 


r (cm- 1 ) 


q 




Concentration 


382.6 


325 


382.6 


325 


2.50 x 10^ u 


8.86 


6.19 


45.4 


-13.2 


4.97 x 10 19 


5.15 


3.66 


46.4 


-13.2 


6.98 x 10 18 


3.50 


2.89 


39.4 


-15.0 



Doped Si (cm 3 ) 


SQ (cm -1 ) 


SQmax 


(cm" 1 ) 


Concentration 


382.6 


325 


382.6 


325 


2.50 x 10^ u 


-3.6 


-1.8 


-3.4 


-2.27 


4.97 x 10 19 


-0.3 


-0.8 


-0.19 


-1.08 


6.98 x 10 18 


0.2 


-0.8 


0.29 


-0.99 
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FIG. 4: Raman lineshapes vs. scattering wavelength for one 
thermally diffused boron doped Si(100) sample with carrier 
concentration 2.50 x 10 20 cm -3 . The experimental data (dis- 
crete points) were then fitted with the modified Fano function 
(solid lines) with a least-squares computer fit which used F 
= 9.39 cm -1 and q as an adjustable parameter. The spectra 
have been shifted vertically with respect to one another. 



B. SOI Samples 

Surface molecular doping of intrinsic silicon-on- 
insulator (SOI) Si(100) samples has been investigated 
and molecular monolayers on hydrogen terminated sil- 
icon for various thicknesses (20 — 200 nm) have been 
formed. Two molecular dopant precursors have been uti- 
lized; allyboronic acid pinacol ester (boronic acid), which 
forms a covalent Si-C bond, and decaborane (B10H14), 
which forms a hydrogen Si-H bond [ll|, [l2[. The sam- 
ples underwent rapid thermal annealing (RTA) for vari- 
ous temperatures and time scales and the carrier concen- 
trations were then calculated by sheet resistance mea- 
surements. RTA is a manufacturing process which heats 
silicon samples to high temperatures (~ 1100°C) on a 
timescale of several seconds. During RTA, the Si-C or 
Si-H bonds are broken and the remaining boron atoms 
become active and diffuse into the surface of the silicon. 

Two SOI samples (75 nm and 25 nm thickness) which 



TABLE III: Parameters obtained by fitting experimental first- 
order Raman spectra with the modified Fano function JT]) and 
holding T = 9.39 cm -1 . Spectra were obtained for a single 
carrier concentration of 2.50 x 10 20 cm -3 . 
Excitation Wavelength (nm) q SQ (cm -1 ) SQ ma x (cm -1 ) 



632.8 


1.47 


-7.8 


-1.41 


532 


3.13 


-7.8 


-4.80 


442 


8.05 


-5.8 


-4.63 


382.6 


27.1 


-3.8 


-3.45 


325 


-23.9 


-1.8 


-2.19 



were doped with the boronic acid precursor and under- 
went RTA (1050° C, 500 s) produced the highest carrier 
concentrations. The calculated carrier concentrations 
were 1.54 x 10 21 cm -3 and 5.20 x 10 21 cm -3 respec- 
tively. The measured intrinsic SOI (75 nm thickness) 
carrier concentration was 8.63 x 10 17 cm -3 . 
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FIG. 5: Raman lineshapes vs. scattering wavelength for three 
molecularly doped SOI Si(100) samples with different carrier 
concentrations. The experimental data (discrete points) were 
then fitted with the modified Fano function (solid lines) with 
a least-squares computer fit which used T and q as adjustable 
parameters. There is no significant change, regardless of the 
carrier concentration. 

The Raman data obtained for the three molecularly 
doped samples are shown in Fig. 5. The experimental 
data (discrete points) were then fitted with the modified 
Fano function (solid lines) with a least-squares computer 
fit which used q and T as adjustable parameters. No sig- 
nificant changes in the Raman lineshape were observed, 
despite the large difference in measured carrier concen- 
tration. 



IV. DISCUSSION 

The parameter q has a strong dependence on the ex- 
citation wavelength as can be seen in Fig. 4. However, 
a transition can be observed in the Raman lineshape be- 
tween the 382.6 nm and 325 nm excitation wavelengths. 
The Raman lineshape for the 325 nm excitation, with 
T held constant, deviates from the theory proposed by 
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Cardona et al. [Ej]. 

One explanation for the divergence is that the car- 
rier concentration is inhomogeneous. The inhomogene- 
ity is not uncommon at such high carrier concentrations, 
which are close to the solubility limit, and would cause 
light with different penetration depths to probe regions 
of different Fermi energies [5j. However, this explanation 
seems unlikely due to the fact that the other two ther- 
mally diffused boron silicon samples, which are at lower 
carrier concentrations, also show this deviation from the- 
ory and a transition between the 382.6 nm and 325 nm ex- 
citation wavelengths. Additionally, all molecular doped 
SOI samples showed no change in the Raman lineshape 
for the 325 nm excitation, despite having the highest car- 
rier concentrations. These results may indicate a tran- 
sition in the silicon band structure from bulk to surface 
bands. 

A second possible explanation for the divergence is to 
consider the penetration depth and electronic excitation 
spectrum for silicon (Fig. 6) 0, BJ, [H, EE E3, • The 
previous work done by Cardona et al. can be consid- 
ered an investigation of the micron region (450 nm — 650 
nm), whereas the present work considers the submicron 
region (325 nm — 632.8 nm). The approximation made 
by Cardona et al. of the combined density of states for 
the continuum of inter- valence band transitions as one- 
dimensional critical points at Er 1 = 3.4 eV and Er 2 = 4.2 
eV leads to a semi-quantitative comparison, since those 
expressions are only valid in the immediate neighborhood 
of the resonant gap [4|, [5[ . 



the Raman scattering. 

By converting the energies of the critical points to exci- 
tation wavelengths, an inflection point in the logarithmic 
plot of the penetration depth in silicon can be seen at 
364 nm. This point is equal to 3.4 eV in energy, which 
corresponds to the first critical T point (EpJ at k = 
in the band structure of silicon. The 382.6 nm (3.24 eV) 
excitation is below the critical point and follows the pro- 
posed theory by Cardona et a/., while the 325 nm (3.8 
eV) excitation is above the critical point and enters a new 
regime of Raman scattering. 
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FIG. 7: Plot of q _1 vs. excitation wavelength for thermally 
diffused boron silicon sample with carrier concentration 2.50 
x 10 20 cm -3 . A power law dependence is observed, which 
deviates from the linear dependence prediction of Cardona et 
al. (noted in outlined box). 
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FIG. 6: Penetration depth vs. wavelength and electronic ex- 
citation spectrum of silicon (inset adapted from [l9j]). An in- 
flection point in the logarithmic plot of the penetration depth 
can be seen at 364 nm. This point is equal to 3.4 eV in energy, 
which corresponds to the first critical F point (Er x ) at k = 
in the band structure of silicon. The 382.6 nm (3.24 eV) 
excitation is below the critical point and follows Cardona's 
theory, while the 325 nm (3.8 eV) excitation is above the 
critical point and enters a new regime of Raman scattering. 

At higher energies, greater than the first critical point 
(E Fl = 3.4 eV), one must consider the anisotropic valence 
bands and calculate the detailed band structure over the 
whole Brillouin zone via a 3D E — k plot to account for 
the inter- valence band transitions that will play a role in 



Additionally, Cardona et al. found a linear dependence 
in the plot of q _1 vs. excitation wavelength for a heavily 
p-type doped silicon sample (1.6 x 10 20 cm -3 ), whereas 
a power law dependence is found for a thermally diffused 
boron silicon sample with carrier concentration 2.50 x 
10 20 cm -3 (Fig. 7) [4]. Cardona et al. probed a smaller 
range of experimental data points, which would appear 
to be linear. However, when analyzing the full spectrum, 
the dependence is clearly non-linear. 



V. CONCLUSIONS 

Based on the experimental and theoretical analysis by 
Raman spectroscopy, for a wide range of excitation wave- 
lengths (325 — 632.8 nm), we have studied the effect of 
inter- valence band transitions on the zone center optical 
phonon and observed Fano interference in heavily p-type 
doped silicon. We were able to fit the asymmetry in 
the one-phonon Raman lineshape by utilizing a modi- 
fied Fano function and have identified a transition above 
and below the one-dimensional critical point (Er x =3.4 
eV) in the electronic excitation spectrum of silicon. By 
analyzing the logarithmic plot of the penetration depth 
in silicon vs. excitation wavelength, we have identified 
an inflection point which corresponds to the first critical 
point Er x and a change in the inter-valence band transi- 
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tions. Energies higher than the first critical point cause 
the system to enter into a new regime of Raman scat- 
tering. As doped silicon devices are scaled down in size, 
Raman spectroscopy is an important tool to study both 
bulk and surface scattering. 
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